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ABSTRACT 



NOAA-7 AVHRR data from April 19S2 was used to perform global-scale analysis 
of aerosol particle characteristics. Mahony's (1991) water vapor correction was incor- 
porated into the AVIIRR multichannel satellite data processing technique used by frost 
(19SS). Channel 4 .'Channel 5 brightness temperature difference was used as an estimate 
of water vapor content in an air column. Greatest measured water vapor content was 
at the equator, decreasing toward the poles. Applying the correction reduced the aerosol 
particle size index, resulting in an increase in aerosol optical depths. In addition, the high 
particle size index over the low latitudes in the southern hemisphere noted by Frost, was 
reduced significantly after applying the correction. A comparison was made between 
Frost's values for the variable two-term Hcnyey-Grecnstcin phase function and those 
derived with the water vapor correction. The basic shape of the curve is similar to that 
determined by Frost. 
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I. INTRODUCTION 



A. PURPOSE 

The availability of satellite-measured radiances provides an effective and relatively 
inexpensive means for determining aerosol particle distribution on a global scale. Satel- 
lites provide extensive coverage of the earth's surface with sufficient temporal periodicity 
to allow for the study of large-scale global effects while at the same time providing 
small-scale local measurements desirable for military applications. Data from 
meteorological satellites require careful processing in order to distinguish the much lower 
radiances associated with aerosol optical depth variations (Durkee 1984) from the 
greater radiances associated with cloud albedos, satellites were designed to measure. 

Improved aerosol optical depth retrieval algorithms have facilitated ongoing efforts 
to establish a global aerosol particle data base. Frost (19SS) developed an algorithm to 
derive aerosol particle characteristics from the Advanced Very High Resolution 
Radiometer (AVHRR) on NOAA polar orbiting satellites. Aerosol particle size index 
can be derived by a comparison of the sensitivities of two adjacent spectral windows and 
so detect aerosol scattering. Frost examined the performance of a variable phase func- 
tion as an improvement to the constant phase function used previously. Mahony (1991) 
examined the effects of water vapor on aerosol optical depth retrieval. An estimate of 
water vapor content in a column of air can be derived by comparing the brightness 
temperature of two AVHRR windows in the thermal infrared region, lie concludes that 
this effect can be compensated for, eliminating an error of 15% in the phase function. 

The purpose of this study is to incorporate Mahony's improvement into the algo- 
rithm used by Frost, which did not account for water vapor appreciably affecting re- 
trieved radiances. The April 19S2 data examined by Frost will be re-evaluated to 
compensate for water vapor. The variable two-term I Ienyey-Greenstein phase function 
recommended by Frost will be examined in light of the new, water vapor corrected par- 
ticle size phase function and compared to the values he generated. 

B. MOTIVATION 

Aerosols consist of a variety of solid or liquid material dispersed in the atmosphere; 
examples include volcanic ash, smoke, and pollen. These particles vary in their distrib- 
ution and concentration in the atmosphere, which can impact local and global radiation 
budgets. The extreme concentrations of smoke particles resulting from the oil fires in 
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Kuwait created a layer that reduced daylight to night in the region. Even modest con- 
centrations have demonstrated an impact in models of the earth's radiation budget 
(Weiss et al. 1974; and Coakley 1976). A Soviet study of longwave radiation in the 
tropical Atlantic ocean (Zaitseva 1976) linked short-term temporal changes in temper- 
ature to changes in aerosol particles. 

Charlson et al. (1987) proposed that the major source of cloud condensation nuclei 
(CCN) over remote unpolluted oceans appears to be a waste product of phytoplankton, 
aqueous dimcthylsulphide gas (DMS). A complex chain of chemical and biological 
processes result in DMS passing through the sea surface. Once in the atmosphere, it is 
oxidized to form sulfate aerosol particles which act as CCN. Varitions in phytoplankton 
populations with the associated changes in CCN particles could alter the earth's heat 
budget. 

Of direct importance to the military is the impact of aerosols on electro-optical de- 
vices. The .4 to 1.0 fim waveband is considered the optimum window for lasers as it 
avoids ozone absorption at shorter wavelengths and water vapor absorption at longer 
wavelengths (Bloembergcn et al. 1987). Aerosol scattering is a significant factor in this 
same waveband. 
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II. THEORY 



A. DIFFUSE RADIANCE 

Knowing that optical depths are small over the ocean and that single scattering is 
dominant (Pfeil 1986), and assuming an optically thin atmosphere. Frost (19SS) used the 
following simplified version of the Radiation Transfer Equation (RTE): 



L = diffuse radiance, 

4> = azimuth angle, 

<y 0 = single scatter albedo, 

F 0 = solar irradiance, 

H = cos 0, ( 0 = satellite zenith angle), 
p = scattering phase function, 

O = single scatter angle, 

<5 = optical depth. 

Griggs (1975) and Durkee et al (1986) demonstrated this linear dependence of radiance 
on optical depth. 

There are three primary sources of satellite-detected radiance; aerosol scatter ( L A ), 
Rayleigh scatter ( L R ), and sea surface reflectance ( L s ). This relationship can be expressed 
as 



. Rayleigh scatter occurs when the size of the aerosol particle is much smaller than the 
wavelength of the incident radiation and, though strongly wavelength dependent, it does 
not vary spatially. Sea surface reflectance is very small for red wavelengths, the albedo 
is 0.5%, and zero for radiation with wavelengths longer than 0.7/un, except where the 
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where, 



L = L a + L r + L s 



2.2 
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sun-earth-satellite geometry results in sunglint (Ramsey 1968). Therefore, by avoiding 
data contaminated by sunglint, total radiance can be approximated by the sum of the 
radiances due to Rayleigh and aerosol scatter (L ^ L A + L R ). By subtracting the compo- 
nent due to Rayleigh scatter, aerosol particle radiance ( L A ) will remain: 

L A = Ui,^4,)-L R ~^^-\mS A . 2.3 

where <5„ is the aerosol particle optical depth. 

With Rayleigh scatter forming a baseline radiance, aerosol scatter is the dominant 
source of upwelled radiance variations ( AL ~ AL a ). Eqn 2.3 indicates satellite measured 
radiance is directly related to the changes in the aerosol particle characteristics of a> 0 , 
p(Q) and d A . It is apparent that a knowledge of p(0) variations is required to determine 
a more accurate <5^ from L A . 

B. OPTICAL DEPTH 

Optical depth is a measure of the scattering and absorption material present in the 
atmosphere. It is determined by vertically integrating the extinction coefficient (<?„,), 
from the earth's surface to the top of the atmosphere (H), expressed by 



<5 = 



’ H 

° ext dz- 

J 0 



2.4 



Extinction is due to both absorption {o cis ) and scattering (o Jccr ) of radiant energy; how- 
ever, in the marine environment, absorption can be neglected, leaving extinction due to 
scattering alone ( a„, » a sc „) such that 



(5 = 



H 



® scat d-- 



2.5 



The amount of scattering is determined by the size of the particles, their scattering 
efficiency, and the size distribution of particles in the column given by 



-r 

J o 



7V' 2 0 



t^scat' 



(m, /, r)n{r)dr, 



2.6 



where. 
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nr 2 = particle cross sectional area, 

Q, ca ,(m. )., r) = scattering efficiency, 
m = complex index of refraction, 

/ = wavelength, 
r = particle radius, 

n(r) = number distribution of particles. 

Aerosol particles are most efficient in scattering radiant energy at wavelengths near 
their radius, becoming less efficient as the difference between the two increases. Any 
changes in the size distribution will produce a corresponding change in <5^ and in the 
observed radiance L A . The variations in radiance between two wavelengths, specifically 
red and near-infrared, can be quantified by the ratio (Durkce ct al. 1991) 

_ {L A )red 

12 (LAir ~ [o> 0 />(O)<5A (> ' 

This ratio (S 12 ) is called the particle size index and is sensitive to variations in the slope 
of the particle size distribution, increasing as the slope increases. Small particles have the 
greatest influence if the slope is steep, while large particles arc most important for more 
gradual slopes. Generally, (L A ) rtd , due to small particle scattering, is greater than {L A ) n , t 
resulting in an S l2 ratio larger than 1 (Mahony 1991). 

C. WATER VAPOR ABSORPTION 

The dominant influence on retrieved radiances in the visible and near-infrared 
wavelengths is aerosol scattering. Though water vapor effects arc not as significant as 
aerosol scatter, they can result in an over-estimation of S 12 which could be interpreted 
as a change in aerosol characteristics or on the other hand, could mask a significant 
change. As the desired result is the estimation of aerosol particle characteristics from 
satellite observations, inclusion of water vapor effects will improve the accuracy of re- 
sults. 

In Figure 1, the effect of water vapor can be seen in the transmittance (r) for 
AVHRR channels 1 and 2. Note that channel I r is greater than channel 2, indicating 
less absorption is taking place. Therefore, L 2 would tend to decrease more due to ab- 
sorption by water vapor than L, leading to the over estimation of S u . 
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Figure 1. Atmospheric transmittance due to outer vapor: Transmittance is 

shown with NO A A & AVI IRK channels 1 and 2 response functions 
(Mahony 1991). 



D. BRIGHTNESS TEMPERATURE DIFFERENCE 

An estimate of the water vapor content in an air column can be made using the 
AV11RR 10 to 13 ft m window region. In channels 4 and 5, at wavelengths of 1 1 and 12 1 
nm respectively, absorption by water vapor is dominant while aerosol scattering is rela- 
th ely insignificant. The benefits of using these channels is the colocation and temporal 
linkage with channels 1 and 2, providing a water vapor snapshot along the same path 
as the aerosol optical depth (<5^) measurements. Although not the ideal sensor windows 
for water vapor content, their use is justified within the accuracy of the approximations 
made (Mahony 1991). 
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The brightness temperature (T) is determined from the radiance received, due to 
direct transfer, from the ocean surface to the satellite sensor. Therefore, the 
transmittance (r) is directly proportional to the retrieved T. Only absorption reduces the 
radiance reaching the sensor; scattering is relatively insignificant. 

Dalu et al. (19S1) proposed a method for determining sea surface temperature using 
brightness temperature difference corrected for air column water vapor content. After 
manipulation of Dalu's formula, the vertically integrated water vapor density (u-) is 
provided by 



w = A(T 4 — T 5 ) cos 6. 2.S 

where A is a constant. Using A = 19600g°A' -1 nr 2 , a scanning angle 0 = 0, and a variety 
of temperature and humidity profiles, Dalu (1986) applied this relationship to a radiative 
transfer model. The resulting line graph (Figure 2) has a correlation coefficient of 
R = 0.99, with an error of ± 1.5 A" m~ 2 in density. 

Mahony determined that the water vapor influence on the particle size index was a 
function of w and S l2 (dry), the latter being the magnitude of S l2 when w = 0. The 
equation for the curve that best fit his data is 



Sn(d r y) = 



Si 2 

1 4- 0.0332b w 



where S l2 and w are determined from measured radiances. 



2.9 



E. VARIABLE THASE FUNCTION 

Based on Frost's work, a variable single scatter phase function was used, having 
been determined to be more accurate in its optical depth estimation. Higher ratios, which 
are due to smaller particles, have more backscatter than the lower ratios, due to larger 
particles. This provides the basic shape necessary for the phase function needed to esti- 
mate optical depth. 

The two-term Henycy-Greenstein scattering phase function (Lenoble 19S5) was 
used. 
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Figure 2. Total water vapor content as a function of brightness temperature differ- 
ence (l)alti 1986). 



/>(©) = 



(1 ~ MQd -g 2 2 ) 

( I + g? - 2g, cos 0) (I + g 2 2 4- g 2 cos 0) 



2.10 



where I V is a weighting factor, g, ami g 2 are asymmetry curve factors, and 0 is the scatter 
angle. All three of these factors are functions of .S u written as 



IF = 1.2 -(0.25 x 5, 2 ), 


2.11 


g, = 1.04 -(0.08333 x S 12 ), 


2.12 1 


g 2 = 1.2 -(0.58333 x S 12 ). 


2.13 



From the equations 2.1 1 through 2.13, the phase function's dependence on S u is appar- 
ent. The marine-type phase function describes large particle distribution, low S 12 ratios; 
while the rural-type phase function follows small particle distribution, high S l2 ratios. 
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As the ratios change, the factors change, resulting in a varying phase fuction from one 
similar to a marine phase function to one similar to a rural phase function. Figure 3 
shows some of the variable phase functions />(©) calculated by Frost (1988) from the 
indicated S 12 values. 
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Figure 3. Aerosol particle si;e index and associated phase functions: (Frost 19S8). 
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III. PROCEDURES 



A. DATA AND PROCESSING 

The data examined was collected by the NO AA 7 polar orbiting satellite whose av- 
erage altitude was 833 km. This sun-synchronous platform crossed the equator at ap- 
proximately 0230 and 1430 local standard time (north- and south-bound, respectively). 
Of the five spectral bands or channels on the AVHRR instrument; channels 1, 2. 4, and 
5 were used in measuring upwelling radiance. Table 1 provides the bandwidths for all the 
channels (Jensen 1986) plus a weighted average of the solar irradiance (F 0 ) over the 
bandwidth used for channels 1 and 2 calculations. 



Table 1. NOAA 7 AVHRR CHANNEL BANDWIDTHS. 



Channel 


Bandwidth 


Solar Irradiance 


i 


0.58-0.68 jam 


154.3536 mWcnr 2 


2 


0.725-1.10 fim 


103.0008 mWcnr 2 


3 


3.55-3.93 f un 




4 


in. 5- 11.3 urn 




5 


1 1.5-12.5 /.un 





The archives of the National Center for Atmospheric Research (NCAR) in Boulder, 
Colorado, provided AVHRR Global Area Coverage (GAC) data for two periods in April 
1982; the 5th - 10th and the 20th - 25th. April 1982 data was analyzed for comparison 
with the work done earlier by Frost (1988). 

All the data was processed on the Interactive Environmental Digital Analysis Lab- 
oratory (IDEA LAB) located at the Naval Postgraduate School, Monterey, CA. The 
IDEA LAB produced enhanced images and graphs from the processed data. 
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B. ANALYSIS 



The data analysis scheme was based on that used by Pfiel (1986) and Frost (1988). 
The water vapor algorithm used by Dalu (1986) was incorporated to account for ab- 
sorption by water vapor, as recommended by Mahony (1991). Also accounted for in the 
algorithm were changes to the solar irradiance due to variations in the sun-earth dis- 
tance. During preprocessing, pixels located at latitudes greater than 70° or those over 
land were discarded. 

To remove sunglint, an empirical method using viewing geometry' was applied. This 
method, a weighted sum of relative zenith and azimuth angles for the sun and satellite, 
identified approximately the same contaminated areas as the Avian sunglint routine, 
which is based on the Cox and Munk (1954) method. 

Ozone absorption effects occur in the red and near infrared wavelengths with the 
greatest absorption at a wavelength of 0.69 nm (Fleagle and Businger, 1980), which is in 
the channel 1 window. This effect was inadvertently compensated for twice by Frost. The 
error has been corrected for this study. 

Each pixel was then examined for clouds as only clear pixels were to be processed. 

Cloud removal was performed by the same method used by Frost: 

1. Test for bad pixels -- If channel 2 or channel 4 radiance values are less than zero, 
after Rayleigh radiance was removed, then the pixel is classified as "bad" and not 
processed. 

2. Test for high clouds -- A pixel was considered to contain high clouds if its channel 
4 brightness temperature was less than 213° K. 

3. Test for low clouds — A pixel was considered to contain low clouds if its channel 
2 albedo was greater than 40%. 

4. Preliminary test for clear pixels -- A channel 1 'channel 2 albedo was calculated and 
if that ratio was greater than 1.5, the pixel was tentatively classified as clear. 

5. Final test for clear pixel — The channel 2 albedo of the pixel tentatively identified 
as clear was compared with the four pixels above, below, to the right and to the 
left. All the pixels must fall with in five counts of each other for the pixel to be 
classified as clear, otherwise it was counted as partlv cloudv. A count is equivalent 
to 0.1069% albedo (Kidwell 1986). 
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The clear pixels were then processed for a variety of parameters. The list is basically 
the same as Frost's, the most significant difference is the inclusion of two parameters 

with the water vapor correction, the terms include the word "dr}’" in their name. 
Channel T Channel 2 ratio, 

Dry Channel 1 Channel 2 ratio, 

Aerosol optical depth, 

Dry Aerosol optical depth, 

Average phase function of Channels 1 and 2, 

Channel 4, Channel 5 brightness temperature differences, 

Number of clear pixels, 

Number of partially cloudy pixels, 

P parameter. 

The P parameter is used to examine the phase function calulated by Frost. 

The phase function p(6, S u ) is an important factor in the determination of aerosol 
optical depth. Any changes made in the phase function algorithm, such as the addition 
of a water vapor correction to S 12 , will produce corresponding changes to the calculated 
optical depths. A comparison of this modified phase function to that determined by 
Frost should provide confirmation that the changes made produce physically sound re- 
sults and add to the knowledge of p(0 , S 12 ). 

A new parameter was defined using the simplified RTF (2.1) that is calculated from 
radiance data and can be used to derive p{6) and S A : 

4/uL 4 

P = —4 ~ p(0, S ]2 )S A . 3.1 

'o 

Using the fact that integrating the phase function over a sphere is unity, the following 
form can be used to determine p(6, S u ): 

Jp(0, S n )da = ^Ap(e,S n )5 A dQ = ^APdQ. = 1, 3.2 

where A is equal to the inverse of the optical depth. It is this equation from which the 
phase function will be calculated for comparison to Frost's curves. 
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All of the results were stored and averaged in 1° by 1° boxes; each data point is one 
box. The number of pixels varied from box to box, depending on the number of partially 
cloudy to cloudy pixels found in each square. 
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IV. RESULTS 



A. INTRODUCTION 

Ongoing efforts to determine aerosol particle characteristics have resulted in a 
steady improvement in the algorithm. Water vapor effects were shown to be an impor- 
tant factor in aerosol optical depth retrieval but have not been previously applied to 
satellite data. In this study the latitudinal variation in water vapor content is initially 
examined to give an overview of where the effects could be the most significant. Particle 
distribution characteristics are examined next, both with and without the water vapor 
correction. The data having the correction will be refered to as "dry". Roth the corrected 
and uncorrectcd particle size index and aerosol optical depth will be examined and 
compared over the Pacific Ocean followed by the latitude profiles of the averaged values. 
Finally, the shape of the P parameter curves arc compared to Frost's phase function 
curves. 

B. COLUMN WATER VAPOR 

Figures 4 and 5 show the average water vapor content in an air column to provide 
a general overview of water vapor effects in the atmosphere for the periods of April 5 - 
10 and 20 - 25, respectively. The general structure is that of a maximum in the vicinity 
of the equator with values then decreasing toward the poles. The drop is more rapid in 
the northern hemisphere than the southern hemisphere, the latter is transitioning from 
its summer season when temperatures arc higher and more water vapor is present. The 
greater year-round air and sea surface temperatures at the equator allow more water 
vapor to exist in a column of air. From this it can be seen that the greatest potential for 
error due to water vapor is in the equatorial region. An examination of the particle size 
index before and after applying the water vapor correction will be made next. 
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LATITUDE PROFILE 




T(4) - T(5) 



Figure 4. Brightness temperature difference for April 5 - 10. 
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LATITUDE PROFILE 




Figure 5. Brightness temperature difference for April 20 - 25. 
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C. rARTICLE SIZE INDEX 



Figures 6 and 7 depict the uncorrected S 12 values across the Pacific Ocean for April 

5 - 10 and 20 - 25, respectively. A broad region of high values extends across the low 
latitudes in the southern hemisphere in both figures, indicating the presense of a greater 
number of smaller particles. Low values indicating larger particles can be seen in Figure 

6 across the northern hemisphere. In Figure 7 there are two noticable maxima; one in 
the vicinity of the New Hebrides islands, the second west of South America. The north 
Atlantic appears to be a region of small particles except along the east coast of the 
United States. 

In Figures 6 and 7, bands cross each figure from north-northwest to south- 
southeast. These are a result of the moving of the sensor's aperature to the left and right 
of the satellite's path to view the side regions. Aerosol scatter is anisotropic, favoring 
forward and backward scatter. A satellite is looking at backscatter when it looks directly 
at the earth. But as the sensor looks to the side, it is no longer viewing backscatter, but 
the smaller scatter associated with smaller scattering angles. This results in a decrease 
in the phase function and a corresponding increase in the optical depth. Future work 
will need to correct for this effect. 

Figures 8 and 9 show a comparison of the particle size index profiles for S l2 and 
S i2 {dry) for early and late April, respectively. The S 12 profiles in both figures show a 
broad region of high values, approximately 2.35, from the equator to 40°S; the values 
then decrease toward the poles more rapidly in the southern hemisphere than in the 
northern. The greatest drop can be seen in the southern hemisphere in Figure 8 with a 
minimum value of 2.0. The minimum in Figure 9 is 2.10. Frost's values were 2.0 at the 
equator decreasing to 1.65 moving southward toward the pole; a second smaller decrease 
occurred toward the north pole, dropping to 1.9. His underestimation was due to a 
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Figure 6. Particle size index over the Pacific, April 5-10 

computational error in the algorithm and over-compensating for ozone absorption 
Next, the corrected S n (dry) will be discussed. 
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Figure 7. Particle size index over the Pacific, April 20 - 25 
D. DRY PARTICLE SIZE INDEX 

Figures 10 and 1 1 depict the S l2 (dry) values across the Pacific Ocean, A comparison 
to the uncorrected S l2 values in I igsures 6 and 7 show that the particle size index values 
decreased significantly after the correction was applied. Open ocean regions are now 
uniform with no significant difference between hemispheres. 
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LATITUDE PROFILE 




Figure 8. Comparison of particle si/e index profiles for April 5 - It): the solid line 

is S u , and the dashed line is 
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rigure 9. Comparison of particle si/e index profiles for A pril 20 - 25: the solid line 

is 5 12 , and the dashed line is 5,j(c/a). 




In early April, Figure 10, distinct high value areas remain in the Gulf of Carpentaria 
and along the coast of China. Late April (Figure 11) shows two areas of higher values 
in the northern Atlantic Ocean and along the west coast of North America. The maxi- 
mum west of South America in Figure 7 is easier to distinguish in Figure 1 1 and has 
lower values. The maxima off the Americas appear to be associated with upwelling re- 
gions that are biologically very productive, which could be a source for the greater 
number of small aerosol particles. 

In Figures 8 and 9, it is easy to see the decrease in size index values between S )2 and 
S n {dry). At the equator S l2 (</rp) is 0.5 less than S, 2 in early April, 0.4 less later in the 
month. From the equator to the poles, the difference between S 12 and S i2 (dry) decreases 
to approximately half the values at the equator in both figures. In both instances the 
overall shape of the S l2 (dry) curve shows a minimum near the equator gradually in- 
creasing as you move toward the poles, which is the reverse of Frost's results where 
values decreased toward the poles. Though a more uniform profile was expected, the 
peaks that exist in the corrected profiles of both figures correspond to the latitudes of 
the maxima depicted in Figures 10 and 1 1. 

E. OPTICAL DEPTH 

The views of the Pacific in Figures 12 and 13, show the southern hemisphere with 
low optical depths and the northern hemisphere with very' high optical depths. The 
largest feature is the wide area of high values, 0.45 to 0.50, that occurs along the 
Mexican Central American west coast extending westward into the Pacific. This feature 
is the result of the El Chicon volcano (17.3°A\ 95.2°I1') erupting on a number of occa- 
sions during the period of 28 March to 6 April. The high values in the Gulf of Mexico 
arc also the result of El Chicon. The movement to the east is possibly due to 
stratospheric winds earning particles ejected to those heights (Frost 1988). In early 



23 




Figure 10. Dry particle size index over the Pacific, April 5 - 10. 

April (Figure 12) the plume extends only halfway across the basin, later in the month 
it appears to approach Japan (Figure 13). 

A second feature can be seen over the East China Sea in Figure 13. This region of 
high optical depth corresponds to a region of large particles in Figure 7. Gobi deseit dust 
blowing off mainland China produced this feature, as noted by Frost (1988). 
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Figure 11. Dry particle size index over the Pacific, April 20 - 25. 

Figures 14 and 15 show a comparison between the optical depth profiles for <5^, and 
d A (dry). The S A profiles both show a maximum value of 0.2 at approximately 20°N, 
which is due to the eruption of the El Chicon volcano. A second peak located in the 
high southern latitudes can also be seen The reason for this maximum is uncertain, 
possibly the result of sea salt from the circumpolar winds around Antarctica. 
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Figure !2. Optica! depth over the Pacific for April 5 - 10. 

F. DRY OPTICAL DEPTH 

Figures 16 and 17 depict d A {dry) across the Pacific Ocean. The high optical depths 
across the northern hemisphere due to El Chicon and Gobi dust have greater areal ex- 
tent and higher values than those in Figures 12 and 13. These features cover regions that 
correspond to the lower S n {dry) regions, indicating larger particles dominate. Figures 14 
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Figure 13. Optical depth over the Pacific for April 20 - 25. 



and 15 are the dry optical depth profiles d A {dry), For the early and late April periods, 
They show that the water vapor correction resulted in an increase in optical depth val- 
ues. Increases were generally on the order of 0.02, but could be as large as 0.06, with 
slightly greater changes occurring at greater distances from the equator. Significant 
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LATITUDE PROFILE 




Optical Depth 



Figure 14. Comparison of optical depth values for April 5 - It): solid lines are S A 

values, dashed are 3J</ry) values. 
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LATITUDE PROFILE 




Optical Depth 



rigure 15. Comparison of optical depth values fur April 20 - 25: solid lines arc <5^ 

values, dashed are values. 
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changes occur in areas where the number of aerosol particles are sufficient to affect the 
optical depth. 

G. PHASE FUNCTION 

Figure 18 shows the P parameter as a function of scatter angle for three different 
particle size indices. A comparison to Frost's plot of the Henyey-Greenstein phase 
function (figure 3) indicates that the basic shape of the P parameter curves are similar 
to that shown by Frost. A smaller particle size index has a larger value in the forward 
direction (between 0° and 60° ) than larger index values. The larger the particle size index 
the closer the curve approaches that of Rayleigh scatter, which is expected as smaller 
particles dominate. Where the scatter angle is between 90° to 140°, the P parameter ap- 
proaches a minimum value for each curve. As the angle increases to 180° the curves 
show an increase in P values. 
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Figure 16. Dry optical depths over the Pacific basin, April 5 - 10. 
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Figure S7. Dry optical depths over the Pacific basin, April 20 - 25. 
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Fig lire 18. I 1 parameter as a function of scatter angle: Values are shown for three 

particle si/e indices; 1.7 (diamonds), 1.9 (triangles), and 2.1 (x's). 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Frost showed that his technique for processing AVHRR data was effective as a tool 
in the analysis of aerosol particle characteristics. Though significant improvements were 
made with his study, sources of error have since been identified requiring futher im- 
provements to his algorithm. The determination that water vapor was a significant 
source of error suggested a re-evaluation of the April 19S2 data in order to assess water 
vapor impact on aerosol optical depth retrieval. 

The results of this study have shown that inclusion of a water vapor correction 
produces a change of approximately 15% in the aerosol particle size index. Without this 
correction, large regions of high particle size index values, which correspond to a 
predominance of relatively smaller particles, were actually regions of lower values 
masked by water vapor effects. This made it difficult to see true areas of small particles 
as they were indistinguishable from the water vapor variations. 

Since the phase function algorithm is a function of S l2 , the modification made for 
water vapor will affect the phase function. It has been shown that the shape of the 
modified phase function is similar to the one determined by Frost. The improved 
parameterization cannot be directly compared to Frost's until the integration (3.2) is 
accomplished. 

The aerosol optical depth has also been shown to be modified by the water vapor 
correction producing an approximate change of 5%. The observed extent and magni- 
tude of areas of high dry optical depths was greater than the optical depths without the 
correction. 

Additional improvements to the algorithm used are still necessary. Multiple scatter 
has not been fully accounted for. The assumption used to allow for single scatter only 
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is valid for optical depths less than about 0.2. This value was exceeded in the regions of 
high optical depth associated with major aerosol particle producing events such as the 
El Chicon eruptions. Errors in the values due to multiple scatter exist and must be cor- 
rected. 

A second area for improvement would be in the determination of S u . Current pro- 
cedures use a ratio of channel 1, channel 2 radiances to calculate the particle size index. 
The preferred method is to use a ratio of channel 1 /channel 2 optical depths, which re- 
quires a better knowledge of the phase function for both channels. It is the optical depth 
which contains information on the aerosol particles and would provide the best estimate 
for S i2 ■ In this study the phase function and optical depths are at the channel 1 wave- 
length. This improvement should eliminate the banding that occurs in the Pacific basin 
photographs which possibly hides features of importance. 

Finally, additional studies of global aerosol particle characteristics are needed. This 
study was limited to a single month. Studies exploring the changes that occur during a 
year and those consisting of a number of years by month or season would increase our 
understanding of the role of aerosol particles in the global climate. 



35 



REFERENCES 



Bloembcrgcn, N'., C. K. Patel, and G. Pake (Chairpersons), 1987: Science and Technol- 
ogy of Directed Energy Weapons. Report of the American Physical Society Study 
Group, Reviews of Modern Physics, 59, S1-S202. 

Charlson, R. J., J. E. Lovelock, M. O. Andrea, and S. G. Warren, 1987: Oceanic 

phytoplankton, atmospheric sulfur, cloud albedo and climate. Nature, 325, 655-661. 

Cox, C. C., and W. Munk, 1954: Measurement of the roughness of the sea surface from 
photographs of the sun's glitter. Journal of the Optical Society of America, 44, 
S3S-S50. 

Coakley, J. A.. 1976: Aerosols and the Earth's radiation budget. Radiation in the At- 
mosphere. IL J. Bollc, ed., Science Press, Princeton, NJ, 472-474. 

Dalu, G., C. Prabhakara, and R. C. Lo, 1981: Improved accuracy of the remote sensing 
of sea surface temperature. In Oceanography from Space, J. F. R. Gower, ed., Ple- 
num, New York, 109pp. 

Dalu, G., 19S6: Satellite remote sensing of atmospheric water vapor. International 
Journal of Remote Sensing, 7, 10S9-1097. 

Durkec, P. A., 1984: The relationship between marine aerosol particles and satellite- 
detected radiance. Ph.D. Thesis, Colorado State University, Fort Collins, CO, US 
ISSN 0067-0340, 124pp. 

Durkee, P. A., D. R. Jensen, E. E. Hindman, and T. II. Vonder Haar, 1986: The re- 
lationship between marine aerosols and satellite detected radiance. Journal of 
Geophysical Research, 91, 4063-4072. 

Durkee, P. A., F. Pfeil, E. Frost, and R. Shema, 1991: Global analysis of aerosol particle 
characteristics. Atmospheric Environment, 25 A- 1 1 , 2457-2471. 

Fleagle. R. W. and J. A. Businger, 1980: An Introduction to Atmospheric Physics, Aca- 
demic Press, New York, NY, 432pp. 

Frost, E. M., 19SS: Global scale estimates of aerosol particle characteristics. Master's 
Thesis, Naval Postgraduate School, Monterey, CA, 1-41. 

Griggs, M.,1975: Measurement of atmospheric aerosol optical thickness over water us- 
ing ERTS-1 data. Journal of Atmospheric Pollution Control Association, 25, 622-626. 

Jensen, J. R.,1986: Introductory Digital Image Processing; A Remote Sensing Approach. 
Prentice-Hall, Englewood ClilTs, NJ, 379pp. 

Kidwell, K. B., 19S6: NOAA Polar Orbiter Data Users Guide. N’OAA N’ESDIS National 
Climatic Data Center, Satellite Data Service Division, 145pp. 



36 



Lenoble, J.,1985: Radiative Transfer in Scattering and Absorbing Atmosphere: Standard 
Computional Procedures., A. Dccpak Pub., Hampton, VA, 300pp. 

Mahony, T. P., 1991: Water vapor influence on satellite-measured aerosol character- 
istics. Master's Thesis, Xaval Postgraduate School, Monterey, CA, 1-38. 

Pfcil, F. R., 1986: Developing a physical basis for an aerosol climatology of the pacific 
ocean. Master's Thesis, Naval Postgraduate School, Monterey, CA, 76pp. 

Ramsev. R. C., 196S: Study of the Remote Measurement of Ocean Color. Final Report, 
TRW, NAS W- 1658. 

Weiss, R., R.J. Charlson, A. P. Waggoner, M. B. Baker, D. Cover, D. Thorsell and S. 
Yuen, 1976: Application of directly measured aerosol radiative properties to climate 
models. Radiation in the Atmosphere IT J. Bolle, ed., Science Press, Princeton, XJ, 
469-471. 

Zaitseva, X. A.. 1976: The spatial-temporal variability in long-wave radiation fields at 
the GATE-74 grounds. Radiation in the Atmosphere . II. J. Bole, ed., Science Press, 
Princeton, XJ, 530-532. 



37 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22304-6145 

2. Library, Code 52 2 

Naval Postgraduate School 

Monterey, CA 93943-5002 

3. Chairman (Code MR'Hy) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey. CA 93943-5000 

4. Professor Phillip A. Durkec (Code MR/De) 2 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5000 

5. Professor Patricia M. Pauley (code MR'Pa) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5000 

6. LT Mary B. Clifford 1 

Fleet Numerical Oceanography Center 

Monterey, CA 93943-5005 

7. Commander 1 

Naval Oceanography Command 

Stennis Space Center 
MS 39529-5000 

S. Commanding Officer 1 

Fleet Numerical Oceanography Center 
Monterey, CA 93943-5005 

9. Commanding Officer 1 

Naval Oceanographic and Atmospheric Research Laboratory 

Stennis Space Center 
MS 39529-5004 

10. Director 1 

Naval Oceanographic and Atmospheric Research Laboratory 

Monterey, CA^93943-5006 



38 



11. Chief of Naval Research 1 

SCO North Quincy Street 

Arlington, VA 22217 

12. Office of Naval Research 1 

Naval Ocean Research and Development Activity 

800 N. Quincy Street 
Arlington, VA 22217 



39 



Thesis 

C5118 Clifford 

cA Effects of water vapor 

and anisotropic scattering 
on aerosol optical depth 
retrieval. 



Thesis 

C5118 Clifford 

c.l Effects of water vapor 

and anisotropic scattering 
on aerosol optical depth 
retrieval. 



